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INTRODUCTION: 


This  study  is  intended  to  investigate  the  physiological  effects  of  antiangiogenic  treatment  at  the 
microcirculatory  level  with  the  intent  to  better  understand  how  such  modulations  of  the  tumor 
environment  can  sensitize  the  tumor  to  conventional  cytotoxic  therapies.  To  this  end,  we  employ  the 
window  chamber  model  to  optically  observe  natural  and  induced  changes  in  tumor  physiology.  A 
number  of  advanced  imaging  modalities  have  been  adapted  to  this  model  in  order  to  facilitate 
longitudinal,  in  vivo  investigation  into  the  parameters  of  interest.  These  include  Doppler  Optical 
Coherence  Tomography  for  the  measurement  of  blood  flow  velocity,  Hyperspectral  Imaging  for  the 
measurement  of  hemoglobin  saturation  and  hematocrit,  and  Confocal  Imaging  for  the  measurement  of 
3-D  vascular  architecture  and  perfusion.  Additionally,  novel  methods  have  been  developed  to 
molecularly  probe  the  tumor  environment.  These  include  Oxygen  Sensitive  Nanoparticles  for  the 
mapping  of  p02  and  the  use  of  Engineered  Cell  Lines  which  express  green  fluorescent  protein  under 
hypoxia.  Together,  these  techniques  will  allow  us  to  better  understand  and  model  the  process  of  tumor 
growth  and  response  to  treatment,  eventually  leading  to  a  better  understanding  of  how  such  changes 
can  be  exploited  to  increase  the  efficacy  of  combinational  therapies. 


BODY: 

The  grant  application's  Statement  of  Work  listed  a  number  of  goals  within  Task  1  to  be  completed  within 
the  first  16  months  of  work.  These  goals  largely  consisted  of  the  development  and  refinement  of 
imaging  modalities  for  the  optical  interrogation  of  tumors  grown  in  a  window  chamber.  While  the 
essence  of  this  project  has  remained  consistent,  we  have  elected  to  employ  the  dorsal  window  chamber 
model  (rather  than  the  mammary  window  chamber  model)  for  the  initial  stages  of  the  project.  While 
the  mammary  window  chamber  model  has  the  advantage  of  allowing  for  orthotopic  tumor 
implantation,  we  felt  that  the  benefits  of  the  dorsal  model  outweighed  this  advantage.  The  ability  to 
use  trans-illumination  to  image  the  dorsal  window  chamber  provides  more  accurate  results  when  using 
the  hyperspectral  imaging  technique,  as  it  does  not  require  complicated  modeling  of  depth-dependant 
tissue  reflectance  and  scattering.  Additionally,  many  of  the  imaging  modalities  would  require  significant 
modification  to  facilitate  the  mammary  model.  For  our  initial  investigation  into  unmodulated  tumor 
progression,  we  felt  that  the  dorsal  model  would  provide  cleaner,  more  easily  processed  data  that 
would  eventually  facilitate  adaptation  to  the  mammary  model,  once  basic  modulations  of  physiological 
parameters  were  thoroughly  modeled  and  understood. 

We  have  also  elected  to  employ  volume  based  optical  measurements  to  a  greater  degree  in  the  initial 
stages  of  assessing  treatment  response.  In  this  design,  tumors  are  grown  in  the  animal's  flank,  and 
physiological  parameters  are  assessed  using  a  fiber  optic  probe.  While  this  technique  is  not  capable  of 
producing  images,  it  has  been  shown  in  earlier  work  that  accurate  values  for  most  of  the  parameters  of 
interest  can  be  non-invasively  serially  measured.  Since  this  technique  is  fast  and  requires  no  surgical 
procedures,  we  can  investigate  time-dependant  treatment  responses  in  a  much  larger  cohort  of  animals 
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than  would  be  possible  using  the  window  chamber  model.  Once  timelines  are  established,  a  window 
chamber  model  will  be  used  in  subgroups  to  collect  more  detailed  and  spatially  quantifiable  data. 

Finally,  the  recent  development  of  oxygen  sensitive  nanoparticles  through  collaboration  with  Cassandra 
Fraser  at  the  University  of  Virginia  has  provided  us  with  an  additional  technique  for  optically  measuring 
extravascular  p02\  We  intend  incorporate  this  invaluable  tool  into  all  of  our  investigations. 

Below  are  presented  the  tasks  originally  proposed  to  be  completed  in  the  first  16  months  of  the  project, 
along  with  discussions  of  accomplishments  or  alterations  to  these  tasks: 

SOW  -  Task  la:  Learn  mammary  window  chamber  surgical  technique  (month  1-2) 

As  discussed  above,  we  have  decided  to  employ  dorsal  window  chamber  and  fiber  probe  based 
techniques,  rather  than  the  mammary  window  chamber  procedure.  I  have  since  become  proficient  at 
both  of  these  techniques  and  have  applied  them  in  a  number  of  preliminary  studies. 

SOW -Task  lb:  Establish  method  for  extracting  fluorescein  concentrations  in  extravascular 
space  using  published  data  by  Dr.  Gregory  Palmer.  Palmer  will  provide  guidance,  (month  1-2) 

This  method  has  been  established  and  will  be  employed  in  future  drug  kinetics  studies. 

SOW  -  Task  lc:  Hemoglobin  saturation  extraction  software  will  need  to  be  adapted  to  mammary 
window  chamber  measurements,  (month  1-2) 


Figure  1,  showing  the  original  hemoglobin  saturation 
reconstruction:  Hemoglobin  saturation  is  represented  as 
grayscale  values,  independent  of  total  hemoglobin  content 
within  the  vessels.  This  representation  has  the  drawback  of 
not  being  able  to  represent  important  variations  in  individual 
vessels'  hemoglobin  (and  thus  oxygen)  delivering  capacity. 


Figure  2,  showing  the  same  tumor,  but  incorporating 
improvements  made  to  the  hyperspectral  processing 
algorithm:  Total  hemoglobin  is  factored  into  the  new 
algorithm  and  is  represented  as  vascular  brightness. 
Hemoglobin  saturation  is  represented  on  a  blue-red  color 
axis,  with  blue  and  red  representing  fully  deoxygenated  and 
oxygenated  blood,  respectively.  The  tumor  appears  in 
white,  while  hypoxia-indicating  GFP  is  represented  by  its 
green  intensity  within  the  tumor. 
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Since  the  dorsal  window  chamber  has  replaced  the  mammary  window  chamber  in  this  aim  of  the  study, 
no  physical  modifications  to  the  hyperspectral  system  were  necessary.  I  have,  however,  made 
improvements  to  the  computational  reconstruction  algorithm  (Figs  1&2),  and  we  are  now  able  to  extract 
total  hemoglobin  values  from  our  data.  Validation  of  this  technique  is  currently  being  performed,  and 
we  expect  to  be  able  to  publish  an  article  describing  this  technique  within  the  next  few  months.  Total 
hemoglobin  will  be  incorporated  into  all  of  our  window  chamber  studies  from  this  point  forward. 

SOW  -  Task  Id:  Practice  imaging  modalities  and  collect  preliminary  data,  (month  3) 

This  task  has  been  completed  and  has  progressed  towards  the  more  detailed  analysis  of  Task  le. 

SOW  -  Task  le:  Collect  HIF-1,  hemoglobin  saturation,  red  blood  cell  flux,  drug  delivery,  and 
vascular  structure  data  for  SA1.  (month  4-15) 


These  parameters  (minus  the 
HIF-1  reporter),  along  with 
redox  ratio,  were 
investigated  in  two  papers 
published  within  the  past 
year2"3.  HIF-1  (as  reported  by 
GFP)  presented  an  obstacle 
to  obtaining  accurate  redox 
ratio  values,  as  the  GFP 
spectrum  overlaps  with  and 
overpowers  the  weak  FAD 
fluorescence  signal  (redox 
ratio  is  measured  as  the  ratio 
of  FAD  to  NADH).  Since  this 
particular  timecourse 

investigated  short  term  physiological  changes  (every  6  hours  for  36  hours),  we  felt  that  redox  ratio  was  a 
better  parameter  to  track,  since  the  long  half-life  of  GFP  might  produce  misleading  quantities  over  such 
a  short  timeframe.  The  data  that  we  collected  yielded  interesting  correlations  among  all  of  these 
parameters  (Fig  3).  Further  data  will  be  collected  and  used  by  a  collaborator  to  develop  detailed  models 
of  vascular  evolution  in  response  to  these  parameters. 

SOW  -  Task  le:  Perform  statistical  analysis  on  data  and  formulate  combinational  treatment 
schedules  for  future  SAs.  These  results  will  be  reviewed  by  a  statistician,  (month  15-16) 

The  project  has  progressed  significantly  ahead  of  schedule  with  regard  to  this  particular  task.  We  have 
already  performed  some  statistical  analysis  on  our  observations,  and  these  results  were  published  in  a 
recent  article  (Fig  4)3.  Furthermore,  we  have  begun  to  develop  advanced  image  processing  algorisms  for 
the  automated  segmentation  and  analysis  of  vascular  architecture.  The  extraction  of  total  hemoglobin 
content  was  one  of  the  first  products  of  this  effort,  and  we  are  well  on  our  way  towards  developing  an 
all  encompassing  algorithm  to  quantify  mean  vessel  density,  tortuosity,  and  mean  vascular  diameter 
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Figure  3,  showing  recently  published  data  regarding  longitudinal  variations  in 
important  physiological  parameters  over  a  36  hour  time  course:  Values  at  each  time 
point  were  averaged  across  all  regions  of  interest.  Hemoglobin  saturation,  flow 
velocity,  vessel  diameter,  redox  ratio,  maximum  velocity,  and  shear  rate  were  all 
shown  to  significantly  change  with  time  (p<0.05). 
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Figure  4,  showing  published  statistical  analysis  of  observed  parameters:  Results 
indicated  a  positive  correlation  between  blood  flow,  hemoglobin  saturation,  and 
metabolic  demand  in  tumors.  But  correlations  are  weak,  suggesting  no  single 
biomarker  can  fully  explain  the  behavior  of  another,  and  further  highlighting  the 
need  for  sophisticated  models  of  tumor  growth  and  adaptation. 


from  our  images.  Such  an 
algorithm  will  allow  us  to 
automatically  extract 
unbiased  values  for  these 
parameters,  circumventing 
the  subjective  and  time- 
consuming  task  of 
performing  manual  image 
analysis. 


KEY  RESEARCH  ACCOMPLISHMENTS: 

•  Developed  appropriate  animal  protocols 

•  Improved  hyperspectral  analysis  to  extract  total  hemoglobin  content 

•  Developed  an  algorithm  for  incorporating  multiple  physiological  parameters  into  a  single  image 
for  better  visual  analysis 

•  Developed  methods  for  incorporating  multiple  imaging  modalities  into  a  single  microscope 
system 

•  Collected  data  over  a  36  hour  time  course  which  showed  significant  fluctuations  in  all 
physiological  parameters  over  that  timeframe 

•  Showed  a  statistically  significant  correlation  among  a  number  of  physiological  parameters, 
suggesting  a  multiparametric  model  of  tumor  progression  is  attainable 


REPORTABLE  OUTCOMES: 

We  have  published  two  papers  over  the  past  two  months  that  demonstrate  significant  progress  towards 
our  goal  of  analyzing  treatment-induced  physiological  changes  and  demonstrating  refined  combinational 
treatment  strategies. 


CONCLUSION: 

The  first  year  of  work  on  this  project  has  produced  highly  optimized  procedures  for  the  analysis  of  key 
physiological  parameters  using  optical  methods.  We  have  further  refined  our  experimental  design  to 
achieve  high  yield  data  which  will  be  incorporated  into  computational  models  for  quick  and  unbiased 
data  extraction.  Additionally,  we  have  begun  to  collect  in  vivo  data  that  has  shown  interesting  results 
over  a  short  time  course.  Statistical  analysis  of  these  data  has  been  published.  We  now  turn  our  focus 
towards  the  final  verification  of  our  methods,  and  continue  on  to  long  time  course  measurements. 
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Computational  methods  for  the  analysis  of  vascular  images  will  be  finalized  in  the  coming  year  and  will 
be  applied  to  the  collected  data.  We  will  soon  investigate  induced  responses  using  cytotoxic  and 
antiangiogenic  therapies  alone.  The  data  gleaned  from  these  studies  will  facilitate  the  exploration  of 
combinational  treatments  in  the  third  year  of  the  project. 
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We  have  combined  hyperspectrai  imaging  with  spectral  domain  optical  coherence  tomography  (SDOCT)  to 
noninvasively  image  changes  in  hemoglobin  saturation,  blood  flow,  microvessel  morphology,  and  sheer  rate 
on  the  vessel  wall  with  tumor  growth.  Changes  in  these  hemodynamic  variables  were  measured  over  24  h  in 
dorsal  skin  fold  window  chamber  tumors.  There  was  a  strong  correlation  between  volumetric  flow  and  he¬ 
moglobin  saturation  (p=0.89,  p  =  9  X  10-6,  N=  15)  and  a  moderate  correlation  between  shear  rate  on  the  ves¬ 
sel  wall  and  hemoglobin  saturation  (p=0.56,  p  =  0.03,  N=  15).  ©  2009  Optical  Society  of  America 
OCIS  codes:  170.4500,  110.4234,  170.2655,  170.3880. 


Hypoxia  is  important  in  regulating  tumor  aggressive¬ 
ness  and  treatment  resistance  and  thus  serves  as  a 
valuable  biomarker  for  tumor  treatment  outcome.  A 
better  understanding  of  tumor  microcirculation 
would  provide  a  basis  for  improved  prognostic  and 
treatment  approaches  for  solid  tumors.  In  tumors, 
formation  of  new  vessels  and  remodeling  and  dropout 
of  existing  vessels  occurs  continuously,  resulting  in 
temporal  variations  in  flow  and  oxygenation  [1].  Dy¬ 
namic  imaging  of  vessel  morphology,  hemoglobin  oxy¬ 
gen  saturation,  blood  flow,  and  shear  rate  on  the  ves¬ 
sel  wall  would  provide  valuable  insight  into  the 
mechanisms  and  distribution  of  cycling  hypoxia  in 
growing  tumors.  For  example,  these  data  will  allow 
for  direct  observation  of  tumor  microvessel  response 
to  changing  hemodynamic  and  metabolic  conditions 
induced  by  cancer  treatments. 

Tumor  hemodynamics  on  the  microvascular  level 
have  traditionally  been  studied  with  invasive  meth¬ 
ods  or  highly  localized  (nonimaging)  methods.  Confo- 
cal  or  multiphoton  microscopy  coupled  with  a  vascu¬ 
lar  contrast  agent  has  been  widely  used  in  animal 
models  to  track  changes  in  vessel  morphology  [2]. 
Blood  flow  in  small  vessels  has  been  limited  to  point 
measurements  using  video-rate  microscopy  and 
single  cell  counting  [1]  or  laser  Doppler  flowmetry 
[3] .  Tumor  oxygenation  has  been  quantified  with  fluo¬ 
rescence  lifetime  imaging  of  a  p02  calibrated  dye  [4] 
or  with  microelectrode  measurements  [3] .  Point  mea¬ 
surements  of  shear  rates  on  vessel  walls  have  also 
been  estimated  with  video  microscopy  [5].  Noninva¬ 
sive  imaging  methods  are  currently  emerging  to  re¬ 
place  these  traditional  techniques.  Speckle  variance 
optical  coherence  tomography  (OCT)  has  recently 
shown  promise  for  monitoring  changes  in  vessel  mor¬ 
phology  with  photodynamic  therapy  in  the  dorsal 
skin  fold  window  chamber  [2] .  Doppler  OCT  imaging 
of  blood  flow  direction  and  velocity  profiles  within 
vessels  [6],  and  hyperspectrai  imaging  of  tumor  he¬ 
moglobin  oxygenation  saturation  [7],  have  also  been 
shown  separately  to  provide  valuable  information  on 
tumor  vessel  function.  However,  a  complete  picture  of 


tumor  hemodynamics  requires  multidimensional  im¬ 
age  data  including  vessel  morphology,  blood  velocity 
profiles  and  direction  of  flow,  and  hemoglobin  oxygen 
saturation. 

We  have  combined  hyperspectrai  imaging  with 
OCT  to  noninvasively  image  changes  in  vascular 
morphology,  hemoglobin  oxygen  saturation,  blood  ve¬ 
locity  and  direction  of  flow,  and  shear  rate  on  the  ves¬ 
sel  wall  with  tumor  growth.  This  combined  micro¬ 
scope  provides  a  wealth  of  information  on  the 
dynamics  of  structural  and  functional  changes  in  tu¬ 
mor  vasculature.  Hyperspectrai  (hemoglobin  oxygen 
saturation)  and  spectral  domain  OCT  (SDOCT) 
(blood  flow  and  vascular  morphology)  arms  were  con¬ 
nected  through  two  separate  baseports  in  an  inverted 
microscope  (Carl  Zeiss  Axiovert  200)  (Fig.  1).  Two- 
dimensional  hyperspectrai  images  were  collected 
with  a  100  W  halogen  lamp  for  transillumination, 
and  detection  was  achieved  with  a  liquid-crystal  tun¬ 
able  filter  (LCTF)  (CRI,  Inc.)  placed  in  front  of  a  DVC 
1412  CCD  camera  (DVC  Company)  [7].  Custom  soft¬ 
ware  was  used  to  tune  the  filter  and  acquire  images 
at  10  nm  increments  between  500  and  610  nm.  Mea¬ 
surements  of  the  dark  offset  and  transmission 
through  a  neutral  density  filter  at  each  wavelength 
were  made  before  each  imaging  session.  Hemoglobin 
saturation  images  were  calculated  by  applying  an  ex¬ 
tension  of  the  Beer-Lambert  law  to  the  wavelength- 
dependent  absorption  at  each  pixel  (assuming  oxy¬ 
genated  and  deoxygenated  hemoglobin  as  the  only 
absorbers),  and  then  solving  for  hemoglobin  satura¬ 
tion  with  linear  least-squares  regression  [7] .  The  sys- 
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Fig.  1.  Schematic  of  the  combined  hyperspectral-SDOCT 
microscope.  L1-L4,  lenses  1-4;  TL,  tube  lens;  BS,  beam 
splitter;  LCTF,  liquid-crystal  tunable  filter. 
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tem,  software,  and  analysis  techniques  have  previ¬ 
ously  been  validated  on  liquid  phantoms  with  an 
accuracy  of  approximately  1%  and  in  vivo  [7] . 

The  common-path  spectral  domain  OCT  arm  con¬ 
sists  of  a  Tiisapphire  laser  source  (Femtolasers)  cen¬ 
tered  at  790  nm  with  a  90  nm  FWHM.  The  interfero- 
gram  is  detected  using  a  custom-made  spectrometer 
with  a  line-scan  CCD  camera  (Atmel,  Aviiva).  The 
OCT  system  is  driven  by  software  that  controls  the 
lateral  scanner  and  performs  data  acquisition  and  ar¬ 
chiving  (Bioptigen,  Inc.).  Three-dimensional  blood 
flow  images  were  collected  with  Doppler  OCT,  which 
measures  phase  changes  due  to  flowing  erythrocytes. 
Phase  changes  are  calculated  from  multiple  A  scans 
collected  at  the  same  position  in  the  volume,  and  can 
be  related  to  flow  velocity  profiles  by  measuring  the 
angle  of  incidence  in  the  three-dimensional  volume. 
Three-dimensional  morphology  was  collected  with 
speckle  variance  OCT,  which  detects  blood  vessels 
from  the  variance  in  the  speckle  pattern  between  se¬ 
quential  B-scan  magnitude  images  [2].  Speckle  vari¬ 
ance  is  advantageous  for  morphology  mapping,  be¬ 
cause  it  is  independent  of  the  angle  between  the 
blood  flow  and  the  incident  beam.  However,  it  is  cur¬ 
rently  not  capable  of  accurately  quantifying  blood  ve¬ 
locity  or  flow  direction. 

Hyperspectral,  Doppler  OCT,  and  speckle  variance 
OCT  images  were  collected  from  4T1  tumors  im¬ 
planted  in  dorsal  skin  fold  window  chambers  in  nude 
mice.  Surgery  and  imaging  were  carried  out  under 
ketamine-xylazine  and  isoflurane  anesthesia,  re¬ 
spectively,  with  the  mice  maintained  at  body  tem¬ 
perature  (conducted  with  institutional  approval  at 
Duke  University).  The  transmission  image  of  the 
window  chamber  morphology  (Fig.  2A)  and  hyper¬ 
spectral  image  of  the  percent  hemoglobin  oxygen¬ 
ation  saturation  in  the  vessels  (Fig.  2C)  were  taken 
with  a  2.5X  objective  (NA=0.12).  Speckle  variance 
OCT  (Fig.  2B)  provides  three-dimensional  vessel 
morphology,  and  Doppler  OCT  (Fig.  2D)  provides 
three-dimensional  vessel  flow  velocities  and  flow  di¬ 
rection.  Speckle  variance  and  Doppler  OCT  images 
were  collected  over  a  1  mm  X  1  mm  area  with  250 
X  125  pixels,  and  1024  pixels  in  the  depth  dimension. 
OCT  images  were  collected  with  a  4x  objective  (NA 
=  0.1)  with  a  1  ms  integration  time  for  each  A  line. 
Eight  repeated  B  scans  were  collected  for  speckle 
variance  OCT,  registered  using  “stackreg”  for  Im- 
ageJ,  and  10  repeated  A  scans  were  colleted  for  Dop¬ 
pler  OCT.  The  flow  profiles  in  Figs.  2E  and  2F  were 
fit  to  a  second-order  polynomial  and  corrected  for  the 
angle  of  incidence  to  provide  velocity  in  millimeters 
per  second.  Maximum  blood  velocity  is  determined 
from  the  peak  of  this  fit,  vessel  diameter  is  deter¬ 
mined  from  the  zero  crossings  of  the  fit,  and  the 
shear  rate  on  the  vessel  wall  is  calculated  from  the 
derivative  of  the  fit  {dvjd r),  assuming  a  Newtonian 
fluid  [8] .  Coregistration  of  speckle  variance  OCT,  hy¬ 
perspectral,  and  Doppler  OCT  images  allow  for  vessel 
morphology,  percent  hemoglobin  oxygen  saturation, 
blood  velocity  and  direction  of  flow,  and  shear  rate  on 
the  vessel  to  be  determined  at  any  point  within  the 
window  chamber. 


Fig.  2.  (Color  online)  Multidimensional  functional  imag¬ 
ing  of  4T1  tumor  vasculature  in  the  dorsal  skin  fold  mouse 
window  chamber.  Transmission  image  of  the  window  cham¬ 
ber  taken  with  the  liquid  crystal  tunable  filter  set  at 
500  nm,  with  a  portion  of  the  tumor  visible  (arrow,  A). 
Scale  bar  is  200  jam.  B,  Speckle  variance  OCT  image  of  ves¬ 
sel  morphology  shown  as  an  en  face  average  intensity  pro¬ 
jection  over  1  mm  depth.  C,  Hyperspectral  image  of  percent 
hemoglobin  oxygen  saturation,  thresholded  for  the 
R-squared  value  of  the  linear  fit  and  total  absorption  value 
at  each  pixel  (background  pixels  are  black).  D,  Doppler 
OCT  image  of  vessel  blood  flow  direction  shown  as  an  en 
face  maximum  intensity  projection  over  depth,  with  red 
vessels  (1)  flowing  toward  the  top  of  the  image  and  blue 
vessels  (2)  flowing  toward  the  bottom  of  the  image.  The 
flow  profiles  at  points  1  (E)  and  2  (F)  were  fit  to  a  second- 
order  polynomial  and  corrected  for  the  angle  of  incidence  to 
provide  velocity  in  millimeters  per  second. 

Changes  in  hemoglobin  saturation,  vessel  morphol¬ 
ogy,  and  blood  flow  with  tumor  growth  were  mea¬ 
sured  in  a  14-day-old  4T1  tumor  in  the  dorsal  skin 
fold  window  chamber  (Fig.  3).  Hyperspectral  and 
Doppler  OCT  data  sets  were  collected  every  6  h  for 
24  h.  The  blood  velocity  and  hemoglobin  saturation 
as  a  function  of  time  are  shown  for  three  vessel  cross 
sections  of  interest  (VOIs)  in  Fig.  3.  The  measured 
hemoglobin  saturation  values  fall  within  the  ex¬ 
pected  range  [9,10].  Note  that  the  p02  of  peripheral 
arteries  in  both  tumors  [9]  and  normal  tissues  [10] 
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Fig.  3.  (Color  online)  Time  series  data  taken  from  a  4T1 
window  chamber  tumor  using  the  combined  hyperspectral- 
OCT  microscope.  The  transmission  image  (A)  is  3  mm 
(height)  X  2  mm  (width)  and  delineates  the  three  VOIs  for 
which  maximum  blood  velocity  (B),  diameter  (C),  and  he¬ 
moglobin  saturation  (D)  are  plotted  (open  circles,  VOI  1; 
open  triangles,  VOI  2;  open  squares,  VOI  3). 
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can  be  substantially  less  than  95%  saturated  owing 
to  loss  of  oxygen  from  the  peripheral  arterial  tree 
(longitudinal  oxygen  gradient).  The  extracted  maxi¬ 
mum  velocities  and  vessel  diameters  are  in  agree¬ 
ment  with  previously  published  values  from  the  win¬ 
dow  chamber  model  [11].  Variations  in  red  blood  cell 
flux  and  p 02  [1],  hemoglobin  oxygen  saturation  [7], 
and  diameter  [9]  have  previously  been  observed  in 
tumors  over  time  scales  ranging  from  one  hour  to 
several  days. 

The  relationship  between  the  maximum  velocity, 
diameter,  flow  rate,  shear  rate,  and  hemoglobin  oxy¬ 
gen  saturation  for  the  same  three  vessels  is  shown  in 
Fig.  4.  Spearman  rank  correlation  analysis  of  the 
three  vessels  over  five  time  points  revealed  a  strong 
correlation  (p=0.94,  p  =  0.018,  N=  5)  between  vessel 
hemoglobin  saturation  and  maximum  vessel  velocity 
for  VOI  2  (Fig.  4B)  and  no  significant  correlation  for 
VOI 1  or  VOI  3  (p  >  0.05,  Figs.  4A  and  4C).  There  was 
no  correlation  between  vessel  diameter  and  maxi¬ 
mum  velocity  or  between  vessel  diameter  and  hemo¬ 
globin  saturation  for  any  of  the  VOIs  (p>0.05).  The 
flow  rate  [Q  =  7t*(Z)/2)2*  Vmax/2(mm3/s)]  was  calcu¬ 
lated  from  the  vessel  diameter  (. D )  and  maximum 
vessel  velocity  (Vmax).  For  all  three  VOIs  and  five 
time  points  grouped  together,  there  was  a  strong  cor¬ 
relation  between  velocity  and  hemoglobin  saturation 
(Fig.  4D,  p=  0.78,  p  =  0.0006,  iV=  15),  and  between  flow 
rate  and  hemoglobin  saturation  (Fig.  4E,  p=  0.89,  p 
=  9XlO-6,  iV=15),  and  a  moderate  correlation  be¬ 
tween  shear  rate  on  the  vessel  wall  and  hemoglobin 
saturation  (Fig.  4F,  p= 0.56,  p  =  0.03,  iV=15).  Shear 
rate  values  are  in  qualitative  agreement  with  previ¬ 
ous  studies  employing  different  techniques  [5] . 

Previous  studies  [1,4]  indicate  that  vessel  remodel¬ 
ing  with  tumor  growth  alters  flow  distributions,  and 


Fig.  4.  (Color  online)  Scatter  plots  of  maximum  velocity 
vs.  hemoglobin  saturation  (O)  and  diameter  vs.  hemoglobin 
saturation  (A)  for  VOI  1  (A),  VOI  2  (B),  and  VOI  3  (C). 
There  was  a  strong  correlation  (B)  between  the  maximum 
velocity  and  hemoglobin  saturation  for  VOI  2  (p= 0.94,  p 
=  0.018,  N=  5),  and  no  correlation  for  any  other  VOI  (p 
>0.05).  However,  when  all  VOI  (O,  VOI  1;  A,  VOI  2;  □, 
VOI  3)  and  time  points  were  grouped  together  (D),  there 
was  a  strong  correlation  between  maximum  velocity  and 
hemoglobin  saturation  (p= 0.78,  p  =  0.0006,  N=15).  There 
was  also  a  strong  correlation  (E)  between  flow  (vessel  cross 
sectional  area  times  \  maximum  velocity)  and  hemoglobin 
saturation  (p=0.89,p  =  9x  10-6,  iV=15)  and  a  moderate  cor¬ 
relation  (F)  between  shear  rate  on  the  vessel  wall  and  he¬ 
moglobin  saturation  (p = 0.56,  p-  0.03,  N- 15)  when  all  VOI 
and  time  points  were  grouped  together. 


the  relationship  between  hemoglobin  oxygen  satura¬ 
tion  and  blood  flow  in  tumors  is  heterogeneous.  How¬ 
ever,  the  techniques  employed  in  previous  studies 
have  a  significant  effect  on  the  results.  The  dual  slit 
technique  suffers  from  inaccurate  estimations  of  the 
length  traveled  by  red  blood  cells,  because  the  length 
of  a  three-dimensional  vessel  is  projected  onto  a  two- 
dimensional  image,  resulting  in  imprecise  velocity 
values.  Red  blood  cell  flux  measurements  are  limited 
to  vessels  smaller  than  —30  pm  in  diameter,  because 
individual  red  blood  cells  are  difficult  to  visualize  in 
larger  vessels.  Also,  previous  studies  measured  intra¬ 
vascular  and  extravascular  p 02,  rather  than  hemo¬ 
globin  saturation.  The  techniques  demonstrated  in 
the  current  study  are  advantageous,  because  flow  ve¬ 
locity  profiles  can  be  accurately  measured  from  a 
wide  range  of  vessel  diameters,  and  morphology, 
blood  oxygen  saturation,  and  blood  velocity  can  be  ac¬ 
curately  determined  without  the  use  of  contrast 
agents. 

We  have  combined  hyperspectral  microscopy  with 
OCT  to  create  a  tool  that  can  dynamically  and  non- 
invasively  monitor  vessel  structure  and  function. 
These  combined  measurements  will  allow  for  unprec¬ 
edented  insight  into  the  relationship  among  vessel 
structure,  hemoglobin  saturation,  blood  flow,  and  hy¬ 
poxic  episodes  in  tumors.  This  research  instrument 
could  ultimately  lead  to  a  more  detailed  understand¬ 
ing  of  the  structural  and  functional  properties  of  tu¬ 
mor  vasculature  and  potentially  become  a  valuable 
tool  for  screening  tumor  therapies  in  animal  models. 
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Abstract.  An  important  feature  of  tumor  hypoxia  is  its  temporal  insta¬ 
bility,  or  "cycling  hypoxia."  The  primary  consequence  of  cycling  hy¬ 
poxia  is  increased  tumor  aggressiveness  and  treatment  resistance  be¬ 
yond  that  of  chronic  hypoxia.  Longitudinal  imaging  of  tumor 
metabolic  demand,  hemoglobin  oxygen  saturation,  and  blood  flow 
would  provide  valuable  insight  into  the  mechanisms  and  distribution 
of  cycling  hypoxia  in  tumors.  Fluorescence  imaging  of  metabolic  de¬ 
mand  via  the  optical  redox  ratio  (fluorescence  intensity  of  FAD/ 
NADFH),  absorption  microscopy  of  hemoglobin  oxygen  saturation, 
and  Doppler  optical  coherence  tomography  of  vessel  morphology  and 
blood  flow  are  combined  to  noninvasively  monitor  changes  in  oxygen 
supply  and  demand  in  the  mouse  dorsal  skin  fold  window  chamber 
tumor  model  (human  squamous  cell  carcinoma)  every  6  h  for  36  h. 
Biomarkers  for  metabolic  demand,  blood  oxygenation,  and  blood 
flow  are  all  found  to  significantly  change  with  time  (p<0.05).  These 
variations  in  oxygen  supply  and  demand  are  superimposed  on  a  sig¬ 
nificant  (p<0.05)  decline  in  metabolic  demand  with  distance  from 
the  nearest  vessel  in  tumors  (this  gradient  was  not  observed  in  normal 
tissues).  Significant  (p<0.05),  but  weak  (r^0.5)  correlations  are 
found  between  the  hemoglobin  oxygen  saturation,  blood  flow,  and 
redox  ratio.  These  results  indicate  that  cycling  hypoxia  depends  on 
both  oxygen  supply  and  demand,  and  that  noninvasive  optical  imag¬ 
ing  could  be  a  valuable  tool  to  study  therapeutic  strategies  to  mitigate 
cycling  hypoxia,  thus  increasing  the  effectiveness  of  radiation  and 

chemotherapy.  ©  201 0  Society  of  Photo-Optical  Instrumentation  Engineers. 
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1  Introduction 

Hypoxic  tumors  are  more  aggressive  than  their  normoxic 
counterparts,  and  are  resistant  to  radiation  and 
chemotherapy,1"3  Thus,  a  major  goal  of  cancer  research  is  to 
gain  a  better  understanding  of  oxygen  supply  and  demand  in 
tumors.  Studies  of  oxygen  supply  and  demand  could  provide  a 
basis  for  improved  prognostic  and  treatment  approaches  for 
solid  tumors.  An  important  feature  of  tumor  hypoxia  is  its 
temporal  instability.  "Cycling  hypoxia”  has  been  described  as 
a  pattern  of  temporal  periodicity  between  hypoxic  and  reoxy¬ 
genated  states.  It  has  been  observed  in  many  tumor  types,4'6 
whereas  studies  measuring  p02  fluctuations  in  normal  tissue 
(muscle)  in  rats  and  mice  have  not  observed  significant 
fluctuations.7,8  One  primary  consequence  of  cycling  hypoxia 
is  upregulation  of  the  transcription  factor  hypoxia-inducible 
factor  1  (HI FI)  activity  to  a  level  that  supersedes  that  of 
chronic  hypoxia.  This  transcription  factor  controls  many  cel- 
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lular  functions  that  exacerbate  treatment  resistance  and  tumor 
aggressiveness.9 

Cycling  hypoxia  could  be  caused  by  changes  in  tumor  vas¬ 
cular  function  (oxygen  delivery)  and  metabolic  demand.  In 
tumors,  formation  of  new  vessels  and  remodeling  and  dropout 
of  existing  vessels  occurs  continuously,  resulting  in  temporal 
variations  in  blood  flow  and  blood  oxygenation.10  Oxygen 
consumption  rate  is  a  highly  dynamic  feature  of  oxygen  trans¬ 
port  in  tumors,  because  small  changes  in  demand  for  oxygen 
create  large  changes  in  the  extent  and  severity  of  hypoxia.11 
Longitudinal  imaging  of  tumor  metabolic  demand,  vessel 
morphology,  hemoglobin  oxygen  saturation,  and  blood  flow 
velocities  would  provide  valuable  insight  into  the  mechanisms 
and  distribution  of  cycling  hypoxia  in  tumors.  For  example, 
these  data  will  allow  for  the  study  of  tumor  response  to 
changing  hemodynamic  and  metabolic  conditions  induced  by 
experimental  and  traditional  cancer  therapies.  Such  informa¬ 
tion  can  also  be  used  to  test  and  refine  theoretical  models  for 
structural  adaptation  in  tumors,  which  have  successfully  pre¬ 
dicted  tumor  response  to  drug  and  radiation  therapy.12 
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Current  methods  for  measuring  oxygen  supply  on  the  mi¬ 
crovessel  level  include  microelectrode  measurements,7  and 
phosphorescence  lifetime  imaging  of  a  p02  calibrated  dye.13 
Microelectrodes  provide  a  sensitive  and  direct  measure  of 
p02,  but  are  inherently  invasive  and  limited  to  point  measure¬ 
ments.  Phosphorescence  lifetime  imaging  provides  informa¬ 
tion  on  the  2-D  distribution  of  oxygen  within  the  tumor,  but 
requires  the  injection  of  a  contrast  agent.  The  use  of  a  contrast 
agent  is  invasive  and  may  require  repeated  exposure  to  the 
dye  for  longitudinal  measurements,  which  further  complicates 
data  acquisition  and  interpretation.  Blood  flow  on  the  mi¬ 
crovessel  level  has  traditionally  been  measured  using  laser 
Doppler  flowmetry,7  the  dual-slit  method,13  or  video 
microscopy.14,15  Laser  Doppler  flowmetry  provides  a  relative 
(not  calibrated)  measurement  of  blood  flow  at  a  single  point, 
and  scanning  laser  Doppler  provides  an  average  flow  over  a 
tissue  area  without  accounting  for  contributions  from  indi¬ 
vidual  vessels.  The  dual-slit  technique  and  video  microscopy 
both  suffer  from  inaccurate  estimations  of  the  length  traveled 
by  red  blood  cells  because  the  length  of  a  3-D  vessel  is  pro¬ 
jected  onto  a  2-D  image,  resulting  in  imprecise  velocity  val¬ 
ues.  Video  microscopy  of  fluorescently  labeled  red  blood  cell 
flux  and  blood  flow  rate  can  be  done,  but  is  only  practical  for 
vessels  less  than  —30  ^m  in  diameter  because  individual  red 
blood  cells  are  difficult  to  visualize  in  larger  vessels  due  to 
absorption  of  fluorescent  light  by  hemoglobin. 

Previous  methods  for  measuring  oxygen  demand  in  tumors 
include  p02  microelectrode  measurements  coupled  with  the¬ 
oretical  simulations  of  oxygen  diffusion,16  an  isolated  tumor 
perfusion  system  combined  with  Fick's  principle 
calculations,1  and  cryospectrophotometric  microtechniques.18 
All  of  these  methods  are  invasive,  which  makes  longitudinal 
in  vivo  monitoring  difficult.  Positron  emission  tomography 
(PET)  of  fluoro-deoxyglucose  (FDG)  has  also  served  as  a 
valuable  marker  for  tumor  metabolism.19,20  Flowever,  PET  im¬ 
aging  is  low  resolution  and  expensive,  and  time  course  mea¬ 
surements  require  repeated  exposure  to  radioisotopes. 

Optical  imaging  of  tumor  oxygen  supply  and  demand  is 
attractive  because  these  methods  are  relatively  low  cost,  high 
resolution,  and  noninvasive,  if  performed  in  a  transparent 
window  chamber  preparation.  Optical  imaging  based  on  en¬ 
dogenous  tissue  contrast  does  not  require  exogenous  contrast 
agents,  and  is  thus  well  suited  for  long-term  monitoring  of 
cycling  hypoxia  in  tumors.  FI yperspectral  imaging  of  hemo¬ 
globin  oxygen  saturation21"23  coupled  with  3-D  Doppler  opti¬ 
cal  coherence  tomography  (OCT)  of  calibrated  microvessel 
blood  flow  rates24"27  allow  for  accurate,  dynamic,  and  nonin¬ 
vasive  imaging  of  oxygen  supply  in  tumors.28  Doppler  OCT 
also  allows  for  visualization  of  vessel  velocity  profiles  at  any 
vessel  cross  section  in  the  tumor,  thus  providing  information 
on  the  shear  rate  on  the  vessel  wall.29  Oxygen  demand  can  be 
optically  monitored  with  the  "redox  ratio,"  which  is  the  ratio 
of  the  inherent  fluorescence  intensity  of  the  metabolic  coen¬ 
zymes  FAD  and  N A D H  (the  primary  electron  acceptor  and 
donor,  respectively,  in  oxidative  phosphorylation.).30  This  op¬ 
tical  redox  ratio  provides  relative  changes  in  the  oxidation- 
reduction  state  in  the  cell.  The  redox  ratio  is  sensitive  to 
changes  in  the  cellular  metabolic  rate  and  vascular  oxygen 
supply.30-34 


The  goal  of  this  study  was  to  dynamically  image  changes 
in  oxygen  supply  and  demand  in  tumors  using  a  novel  com¬ 
bination  of  redox  imaging  (metabolic  demand)  and  hyper- 
spectral  imaging,  combined  with  Doppler  OCT  (oxygen  sup¬ 
ply).  The  results  reported  here  are  the  first  demonstration  of 
combined  noninvasive  imaging  of  oxygen  supply  and  meta¬ 
bolic  demand  in  tumors  in  vivo.  These  data  could  provide 
insight  into  the  sources  of  cycling  hypoxia  in  tumors,  and 
serve  as  the  basis  for  further  studies  that  will  incorporate 
therapeutic  strategies  to  mitigate  hypoxia  in  tumors  for  im¬ 
proved  treatment  response. 


2  Materials  and  Methods 

2.1  Multifunctional  Microscope 

Fluorescence,  hyperspectral,  and  OCT  imaging  were  con¬ 
ducted  on  the  same  microscope  to  noninvasively  monitor 
changes  in  metabolic  demand  and  hemodynamics  with  tumor 
growth.  All  images  were  collected  with  a  4x  objective  (NA 
=0.1).  FI  yperspectral  (hemoglobin  oxygen  saturation)  and 
SDOCT  (blood  flow  and  vascular  morphology)  arms  were 
connected  through  two  separate  base  ports  in  an  inverted  mi¬ 
croscope  (Carl  Zeiss  Axiovert  200).  2-D  hyperspectral  images 
were  collected  with  a  100-W  halogen  lamp  for  transillumina¬ 
tion,  and  detection  was  achieved  with  a  10-nm  bandwidth 
liquid  crystal  tunable  filter  (LCTF)  (CRI  Incorporated, 
Woburn,  Massachusetts)  placed  in  front  of  a  DVC  1412  CCD 
camera  (DVC  Company,  Austria,  Texas).23  Custom  software 
was  used  to  tune  the  filter  and  acquire  images  at  10-nm  in¬ 
crements  between  500  and  620  nm.  M  easurements  of  the  dark 
offset  and  transmission  through  a  neutral  density  filter  at  each 
wavelength  were  made  before  each  imaging  session.  Flemo- 
globin  saturation  images  were  calculated  by  applying  an  ex¬ 
tension  of  the  Beer-Lambert  law  to  the  wavelength-dependent 
attenuation  at  each  pixel  (assuming  oxygenated  and  deoxy- 
genated  hemoglobin  as  the  primary  absorbers,  along  with  a 
tissue  scattering  parameter),  and  then  solving  for  hemoglobin 
saturation  with  linear  least-squares  regression.23  The  system, 
software,  and  analysis  techniques  have  previously  been  vali¬ 
dated  on  liquid  phantoms  with  an  accuracy  of  approximately 
1%  and  in  vivoP 

2-D  fluorescence  imaging  was  conducted  in  the  epi- 
illumination  geometry  with  a  100-W  mercury  lamp.  N A D H 
images  were  collected  with  a  360-nm  bandpass  excitation  fil¬ 
ter  (40-nm  bandwidth)  and  a  390-nm  longpass  dichroic 
beamsplitter.  FAD  images  were  collected  with  a  470-nm 
bandpass  excitation  filter  (40-nm  bandwidth)  and  a  510-nm 
longpass  dichroic  beamsplitter.  Fluorescence  emission  was 
collected  with  the  same  LCTF  and  charge-coupled  device 
(CCD)  described  before  for  the  hyperspectral  arm.  The  LCTF 
was  set  to  480-nm  and  530-nm  emission  for  N  A  D  H  and  FAD 
imaging,  respectively.  Redox  images  were  calculated  by  di¬ 
viding  the  FAD  image  point  by  point  by  the  N A D H 
image.30,34  To  account  for  interanimal  variability  in  the  time- 
course  experiments,  the  redox  ratio  for  each  animal  was  nor¬ 
malized  the  first  time  point  of  the  first  region  of  interest 
within  each  animal.  Note  that  the  redox  image  represents  the 
redox  ratio  of  a  superficial  (—200  to  300  pm  thick)  tissue 
layer  probed  by  the  excitation-emission  wavelengths,35  and 
the  hemoglobin  saturation  images  show  vessels  from  a  simi- 
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larly  superficial  tissue  layer.22,23  OCT  volumes  indicate  that 
this  layer  contains  most  of  the  blood  vessels  in  the  field  of 
view.  The  spatial  resolution  of  hyperspectral  and  redox  imag¬ 
ing  is  defined  by  the  NA  of  the  objective  (5  [xm ). 

The  spectral  domain  OCT  arm  was  illuminated  by  a 
Ti:sapphire  laser  source  (Femtolasers,  Vienna,  Austria)  cen¬ 
tered  at  790  nm  with  a  90-nm  full  width  at  half-maximum. 
The  OCT  interferometer  was  common  path,  utilizing  the  top 
surface  of  the  window  chamber  cover  glass  as  the  reference 
reflection,  and  the  spectral  interferogram  was  detected  using  a 
custom-made  spectrometer  (spectral  resolution  and  coverage 
0.08  and  164  nm,  respectively)  with  a  2048  element  line  scan 
CCD  camera  (Atmel,  Aviiva,  San  Jose,  California).  The  OCT 
system  was  driven  by  software  that  controls  the  lateral  scan¬ 
ner  and  performs  data  acquisition  and  archiving  (Bioptigen, 
Incorporated,  Research  Triangle  park,  North  Carolina).  3-D 
blood  flow  images  were  collected  with  Doppler  OCT,  which 
measures  phase  changes  due  to  flowing  erythrocytes.24  Phase 
changes  were  calculated  from  multiple  A-scans  collected  at 
the  same  position  in  the  volume,  and  were  related  to  flow 
velocity  profiles  by  measuring  the  angle  of  incidence  in  the 
3-D  volume.  This  method  has  previously  been  validated  in 
flow  phantoms.36  Doppler  OCT  volumes  were  collected  over 
a  2x1.5  mm  area  with  250x  125  pixels  (with  a  resulting 
resolution  of  8  and  12  ^m  in  the  lateral  dimensions),  and 
from  2  mm  in  the  depth  dimension  with  1024  pixels,  with  a 
1-ms  integration  time  for  each  A-line.  Ten  repeated  A-scans  at 
each  location  were  collected  at  each  x-y  position  for  Doppler 
shift  estimation.  Cross  sectional  vessel  flow  profiles  were  fit 
to  a  second-order  polynomial  and  corrected  for  the  angle  of 
incidence  to  provide  velocity  in  millimeters  per  second.  Maxi¬ 
mum  blood  velocity  (Vmax)  was  determined  from  the  peak  of 
this  fit,  vessel  diameter  ( D )  was  determined  from  the  zero- 
crossings  of  the  fit,  and  the  shear  rate  on  the  vessel  wall  was 
determined  from  the  derivative  of  the  fit  (dvjdr)  evaluated  at 
the  vessel  wall,  assuming  a  Newtonian  fluid.29  Note  that  in 
vessels  with  diameters  30  to  60  p. m ,  the  velocity  profile  is 
expected  to  be  blunter  than  a  parabola  (this  reflects  the  non¬ 
continuum  behavior  of  blood  rather  than  non-Newtonian 
behavior).25,37-39  Thus,  the  estimated  wall  shear  rate  may  be 
an  underestimate  in  these  calculations.  However,  our  fit  re¬ 
sults  suggest  that  the  departure  from  a  parabola  is  small,  so 
the  second-order  polynomial  fit  is  a  reasonable  approxima¬ 
tion.  The  flow  in  the  vessel  is  the  average  velocity  in  the 
vessel  times  the  cross  sectional  area  of  the  vessel 
[Vmax/2*77*(£>/2)2]. 

2.2  In  Vivo  Imaging 

All  in  vivo  experiments  were  conducted  under  a  protocol  ap¬ 
proved  by  the  Duke  University  Institutional  Animal  Care  and 
Use  Committee.  A  titanium  window  chamber  was  surgically 
implanted  under  anesthesia  (ketamine  100  mg/kg  and  xyla- 
zine  10  mg/kg  IP)  on  the  back  of  athymic  nude  mice  (nu/nu, 
NCI,  Frederick,  Maryland).  A  window  chamber  tumor  was 
established  during  chamber  implantation  by  injecting  20  [xL 
of  a  single  cell  suspension  (8xl06  cells)  of  FaDu  (human 
squamous  cell  carcinoma)  cells  into  the  dorsal  skin  fold  prior 
to  placing  a  12-mm-diam  number  2  round  glass  coverslip 
(Erie  Scientific,  Portsmouth,  New  Hampshire)  over  the  ex¬ 
posed  tissue.  Animals  were  housed  in  an  on-site  housing  fa¬ 


cility  with  ad  libitum  access  to  food  and  water  and  standard 
12-h  light/dark  cycles. 

Redox  and  hemoglobin  saturation  imaging  were  validated 
on  mice  with  a  tumor-free  window  chamber.  Perturbation  ex¬ 
periments  were  performed  with  the  mice  under  ketamine/  xy- 
lazine  anesthesia  (maintained  at  body  temperature)  while 
breathing  air,  then  100%  oxygen  through  a  nose  cone,  and 
then  100%  nitrogen  through  a  nose  cone  at  1-L/min  flow 
rate. 

For  longitudinal  imaging  of  tumor  growth,  animals  were 
anesthetized  with  Isoflurane  (1.5%  with  oxygen)  and  main¬ 
tained  at  body  temperature.  Imaging  began  when  tumors 
reached  1  to  2  mm  in  diameter.  A  Doppler  OCT  volume 
(5-min  acquisition),  followed  by  hemoglobin  saturation  (30-s 
acquisition)  and  redox  (15-sec  acquisition)  imaging  were  col¬ 
lected  every  6  h  for  36  h.  Vessel  morphology  maps  were  used 
to  achieve  spatial  registration  between  each  imaging  session. 
Transition  between  microscope  modes  was  approximately 
1  min.  M  ice  were  allowed  to  recover  after  each  imaging  ses¬ 
sion. 

Note  that  no  corrections  were  made  to  the  data  to  account 
for  the  cardiac  cycle.  All  vessels  of  interest  were  microvessels 
(diameter  <70  /urn)  that  do  not  exhibit  large  changes  in  flow 
due  to  the  cardiac  cycle.  Previous  studies  have  used  Fourier 
analysis  of  blood  flow  and  p02  in  rats  to  show  low-frequency 
fluctuations  (<2  cycles/min)  in  both  tumor  and  muscle  tis¬ 
sues,  which  is  much  lower  than  the  measured  respiration  rate 
and  heart  rate  (~50  to  60  and  300  to  400  cycles/min, 
respectively.7 

2.3  Statistical  Analysis 

Doppler  OCT  and  hyperspectral  data  were  collected  over  the 
36-h  time  course  from  a  total  of  five  mice,  and  redox  images 
were  also  collected  over  the  time  course  from  three  of  these 
mice.  Three  to  six  regions  of  interest  were  chosen  within  each 
tumor  based  on  visual  inspection  of  the  hemoglobin  saturation 
time  course  images.  Vessel  regions  that  demonstrated  fluctua¬ 
tions  in  hemoglobin  saturation  over  the  time  course,  and  that 
were  spatially  distributed  across  the  surface  of  the  tumor, 
were  chosen  as  regions  of  interest.  The  mean  and  standard 
deviation  of  the  hemoglobin  saturation  were  measured  at 
these  vessel  regions  of  interest,  over  an  area  within  the  inner 
diameter  of  the  vessel  and  a  slab  thickness  of  no  more  than 
20-  ytim.  The  Doppler  velocity  profile  was  measured  at  the 
same  vessel  cross  section  of  interest.  The  mean  and  standard 
deviation  of  the  redox  ratio  was  measured  from  a 
25  to  60-^t.m-diam  circular  area  adjacent  to  the  vessel.  The 
redox  area  of  interest  was  chosen  to  exclude  blood  vessels  and 
to  remain  less  than  100  /x m  in  diameter  (the  length  of  oxygen 
diffusion  in  tissue.40  A  total  of  21  regions  of  interest  were 
analyzed  over  the  time  course,  with  14  regions  of  interest  that 
include  redox  ratio  measurements.  One  of  the  animals  for 
which  redox  images  were  collected  included  only  six  time 
points  for  the  four  regions  of  interest  within  that  animal,  due 
to  tumor  hemorrhage  at  the  last  time  point. 

Pearson's  correlation  coefficients  were  used  to  determine 
linear  correlations  between  blood  oxygenation  (hemoglobin 
saturation),  blood  flow  (maximum  velocity,  flow,  vessel  inner 
diameter,  and  shear  rate),  and  metabolic  demand  (redox  ratio). 
Significant  correlations  (/? < 0.05)  are  reported  for  the  pooled 
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data  (all  regions  of  interest  and  time  points  grouped  together), 
for  which  there  were  n  = 143  samples  for  the  blood  oxygen¬ 
ation  and  blood  flow  comparisons,  and  n= 94  samples  for  the 
metabolic  demand  comparisons.  As  described  before,  the  re¬ 
dox  ratio  was  normalized  to  the  first  time  point  of  the  first 
region  of  interest  within  each  animal.  Thus,  blood  oxygen¬ 
ation  and  blood  flow  biomarkers  were  also  normalized  to  the 
first  time  point  of  the  first  region  of  interest  within  each  ani¬ 
mal  for  comparisons  with  the  redox  ratio  only. 

Next,  a  linear  mixed  model  was  used  to  analyze  repeated 
measures.  This  approach  separates  the  interanimal  variation 
from  the  total  variation  (thus  increasing  statistical  power)  to 
determine  whether  the  biomarkers  (i.e„  hemoglobin  satura¬ 
tion,  redox  ratio,  maximum  velocity,  vessel  inner  diameter, 
flow,  and  shear  rate)  change  over  time.  Let  Yijk  denote  the 
value  of  a  biomarker  from  the  ;'th  animal  in  the / th  region  of 
interest  (ROI)  at  the  &’th  time  point.  Assume  Yijk  follows  the 
linear  mixed  model  Yijk=/j.+tk+bij+eijk  (model  1);  where 
fi+tk  is  the  time  effect,  is  the  animal  specific  random 
effect,  and  sijk  are  intra-animal  errors.  The  time  effect  was 
tested  by  the  F-statistic.  A  second  mixed  model  was  used  to 
address  the  association  of  two  variables  adjusting  for  the  co¬ 
variate  time.  Model  2  ( Yljk = +Xljk+tk Ybij Jrsl]k)  includes 
the  covariate  X,  and  the  linear  association  of  biomarker  X  on 
the  behavior  of  another  biomarker  Y  was  tested  with  the 
F-statistic.  Associations  were  tested  between  blood  oxygen¬ 
ation  (hemoglobin  saturation),  blood  flow  (maximum  velocity, 
flow,  vessel  inner  diameter,  and  shear  rate),  and  metabolic 
demand  (redox  ratio)  biomarkers. 

3  Results 

The  results  of  perturbation  experiments  in  a  normal  mouse 
window  chamber  (no  tumor  cells  were  injected)  are  shown  in 
Fig.  1.  To  segment  the  vessels  from  background,  the  hemo¬ 
globin  saturation  images  were  thresholded  by  the  total  absorp¬ 
tion  image.  The  same  mask  was  negatively  applied  to  the 
redox  images.  Note  that  there  are  two  distinct  peaks  in  the 
hemoglobin  saturation  histogram  for  air  and  oxygen  breathing 
[Fig.  1(b)],  due  to  arteriolar  and  venous  blood  oxygenation. 
Nitrogen  and  oxygen  breathing  resulted  in  a  decrease  and  in¬ 
crease,  respectively,  in  both  the  redox  ratio  and  the  hemoglo¬ 
bin  saturation  compared  to  air  breathing,  as  expected.  The 
results  in  Fig.  1  were  confirmed  in  a  second  mouse,  and  are 
consistent  with  previous  reports.23,33  These  perturbation  ex¬ 
periments  indicate  that  NADH  and  FAD  signals  have  been 
properly  isolated  for  redox  imaging. 

The  hemoglobin  saturation  [Fig.  1(a)]  and  redox  [Fig. 
1(b)]  images  collected  from  one  mouse  over  the  36-h  time 
course  are  shown  in  Fig.  2.  Longitudinal  oxygen  gradients  are 
apparent  in  the  hemoglobin  saturation  images,  with  decreased 
hemoglobin  saturation  levels  within  the  tumor  [Fig.  2(a),  at 
the  1]  compared  with  the  tumor  periphery  [Fig.  2(a),  at  the  2], 
These  longitudinal  gradients  have  been  previously  observed 
in  studies  employing  phosphorescence  lifetime  imaging  of 
/?02.41,42  The  redox  images  indicate  that  the  redox  ratio  within 
the  tumor  is  higher  in  areas  with  large-diameter  oxygenated 
vessels  [Fig.  2(b),  at  the  3],  Both  the  hemoglobin  saturation 
and  the  redox  ratio  appear  to  vary  with  time. 

Figure  3  shows  an  en-face  view  of  the  3-D  Doppler  OCT 
volume  (manually  segmented  in  Amira  software,  Mercury 


(a) 


Fig.  1  Hemoglobin  saturation  (HbSat)  and  redox  images  taken  from  a 
nontumor-bearing  window  chamber  of  a  mouse  breathing  oxygen,  air, 
and  nitrogen  (a)  at  1  L/min  through  a  nose  cone.  To  segment  the 
vessels  from  background,  the  hemoglobin  saturation  images  were 
thresholded  by  the  total  absorption  image.  The  same  mask  was  nega¬ 
tively  applied  to  the  redox  images.  Histograms  of  (b)  the  hemoglobin 
saturation  and  (c)  redox  images  show  that  the  hemoglobin  saturation 
and  redox  ratio  increase  with  oxygen  breathing  and  decrease  with 
nitrogen  breathing,  as  expected.  Image  sizes  are  2X  1.5  mm. 


Systems,  Chelmsford,  M  assachusetts)  from  the  zero-hour  time 
point  [Fig.  3(a)]  of  the  same  animal  shown  in  Fig.  2.  Plots  of 
the  cross  sectional  vessel  velocity  profile  from  one  vessel  [ar¬ 
row,  Fig.  3(a)]  are  shown  along  with  the  second-order  poly¬ 
nomial  fit  and  the  R-squared  value  of  the  fit  for  each  time 
point  [Figs.  3(b)- 3(h)].  The  shape  and  peak  of  the  velocity 
profile  changes  with  time,  indicating  changes  in  vessel  diam¬ 
eter,  shear  rate,  and  maximum  vessel  velocity  over  the  36-h 
time  course. 

Figure  4  shows  quantitative  measures  of  the  redox  ratio 
[Fig.  4(b)],  hemoglobin  saturation  [Fig.  4(c)],  vessel  maxi¬ 
mum  velocity  [Fig.  4(d)],  vessel  inner  diameter  [Fig.  4(e)], 
flow  [Fig.  4(f)],  and  shear  rate  [Fig.  4(g)]  over  the  36-h  time 
course  in  six  regions  of  interest  [Fig.  4(a)]  from  the  same 
animal  shown  in  Figs.  2  and  3.  These  plots  indicate  that  all 
variables  change  with  time  within  the  same  region  of  interest 
within  the  same  animal.  All  regions  of  interest  appear  to  vary 
together  with  time  for  the  redox  ratio  [Fig.  4(b)]  and  for  the 
hemoglobin  saturation  [Fig.  4(c)].  However,  each  region  of 
interest  appears  to  vary  differently  with  time  compared  with 
the  other  five  regions  of  interest  for  the  blood  flow  parameters 
[Figs.  4(d)- 4(g)], 

The  spatial  distribution  of  the  redox  ratio  with  respect  to 
the  nearest  vessel  is  shown  in  Fig.  5  for  tumors  and  normal 
tissues.  There  is  a  significant  (/? < 0.05)  decrease  in  the  redox 
ratio  for  points  greater  than  30  /xm  from  the  nearest  vessel, 
compared  to  the  point  closest  to  the  nearest  vessel  (point 
zero).  This  is  in  agreement  with  previous  reports  that  mea¬ 
sured  p02  as  a  function  of  distance  from  the  nearest 
vessel.15,16  However,  the  normal  tissues  showed  no  significant 
increase  or  decrease  in  the  redox  ratio  as  a  function  of  dis¬ 
tance  from  the  nearest  vessel. 
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Fig.  2  (a)  Hemoglobin  saturation  and  (b)  redox  images  taken  from 
one  mouse  over  the  36-h  time  course.  The  time  of  imaging  in  hours  is 
indicated  in  the  lower  left  corner  of  each  hemoglobin  saturation  im¬ 
age.  For  the  purposes  of  this  figure  only,  the  hemoglobin  saturation 
images  were  thresholded  by  the  total  absorption  image,  and  then  by  a 
1 5-/xm-diam  circular  kernel.  The  kernel  excluded  pixels  with  an  ab¬ 
sorption  value  smaller  than  the  average  of  the  15-/xm-diam  kernel. 
The  mask  was  applied  to  hemoglobin  saturation  images  to  isolate 
vessels  from  the  nonvascular  background,  and  the  same  vessel  mask 
was  negatively  applied  to  the  redox  images.  The  tumor  boundary  is 
indicated  by  a  white  dashed  line  in  (b).  Image  sizes  are  2  X  1 .5  mm. 


The  relationship  between  the  biomarkers  for  blood  oxy¬ 
genation  (hemoglobin  saturation),  blood  flow  (maximum  ve¬ 
locity,  flow,  vessel  inner  diameter,  and  shear  rate),  and  meta¬ 
bolic  demand  (redox  ratio)  were  tested  with  the  Pearson's 
correlation  coefficient  (Table  1)  for  the  pooled  data  (all  re¬ 
gions  of  interest  and  time  points  from  all  animals  grouped 
together).  The  correlation  coefficients  (r)  in  Table  1  are 
shown  for  biomarker  relationships  for  which  a  significant  (p 
<0.05)  linear  correlation  was  found.  The  hemoglobin  satura¬ 
tion  is  positively  correlated  with  the  maximum  velocity,  flow, 
and  shear  rate  (HbSat  line  of  Table  1),  and  the  redox  ratio  is 
positively  correlated  with  the  hemoglobin  saturation,  maxi¬ 
mum  velocity,  flow,  and  diameter.  All  statistically  significant 
correlations  are  weak  (r«0.5),  which  indicates  that  many 
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Fig.  3  An  en-face  view  of  the  3-D  Doppler  OCT  volume  (manually 
segmented  in  Amira  software,  Mercury  Systems)  from  (a)  the  zero- 
hour  time  point  of  the  same  animal  shown  in  Fig.  2  (2Xl.5-mm 
imaged  area,  0.6  mm  depth).  Plots  of  the  cross  sectional  vessel  veloc¬ 
ity  profile  from  one  vessel  [arrow  in  (a)]  are  shown  along  with  the 
second-order  polynomial  fit  and  the  R-squared  value  of  the  fit  for  each 
time  point  (b)  through  (h).  The  time  of  imaging  in  hours  is  indicated  in 
the  upper  right  corner  of  each  plot. 


factors  contribute  to  the  variability  of  each  biomarker,  as 
expected.11  The  linear  mixed  model  determined  that  all  biom¬ 
arkers  significantly  change  (p<0.05)  with  time  (Fig.  6),  and 
that  the  vessel  diameter  is  associated  with  the  redox  ratio  of 
the  adjacent  tissue  (p < 0.05)  after  adjusting  for  the  effect  of 
time.  No  other  biomarker  pairs  were  found  to  have  significant 
associations  (p >  0.05)  when  adjusting  for  the  effect  of  time. 

4  Discussion 

Multifunctional  optical  imaging  of  oxygen  supply  and  de¬ 
mand  has  several  advantages  over  traditional  techniques.  The 
combination  of  OCT,  absorption,  and  fluorescence  imaging 
provides  quantitative  information  on  blood  flow,  blood  oxy¬ 
genation,  and  metabolic  demand,  respectively,  in  living  tis¬ 
sues.  This  combination  of  optical  imaging  modalities  is  non- 
invasive,  relatively  inexpensive,  and  fast.  These  attributes  are 
well  suited  for  long-term  repeated  measures  in  vivo.  We  have 
demonstrated  the  attributes  of  this  imaging  platform  for  the 
study  of  cycling  hypoxia  in  tumors  in  the  dorsal  skin  fold 
mouse  window  chamber.  This  is  the  first  demonstration  of 
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Fig.  4  Quantitative  measures  of  (b)  the  redox  ratio,  (c)  hemoglobin 
saturation,  (d)  vessel  maximum  velocity,  (e)  vessel  inner  diameter,  (f) 
flow,  and  (g)  shear  rate  over  the  36-h  time  course  in  (a)  six  regions  of 
interest  from  the  same  animal  shown  in  Figs.  2  and  3.  The  mean  and 
standard  error  are  plotted  for  (b)  the  redox  ratio  and  (c)  hemoglobin 
saturation,  and  the  mean  (d)  through  (g)  is  plotted  for  the  remaining 
variables. 


combined  noninvasive  monitoring  of  oxygen  supply  and 
metabolic  demand  in  tumors  in  vivo. 

All  biomarkers  (hemoglobin  saturation,  redox  ratio,  maxi¬ 
mum  velocity,  vessel  inner  diameter,  flow,  and  shear  rate) 
significantly  changed  with  time  (/? < 0.05)  over  the  36-h  im¬ 
aging  period,  which  supports  previous  reports  that  blood  oxy¬ 
genation,  blood  flow,  and  metabolic  demand  all  contribute  to 
cycling  hypoxia  in  tumors.3  These  variations  in  oxygen  sup¬ 
ply  and  demand  are  superimposed  on  a  decline  in  metabolic 
demand  with  distance  from  the  nearest  vessel  in  tumors  (Fig. 
5),  in  agreement  with  previous  interstitial  p02  microelectrode 
measurements.15,16  Our  results  indicate  a  positive  correlation 
between  the  blood  flow,  blood  oxygenation,  and  metabolic 
demand  in  tumors  (Table  1),  but  the  correlation  between  these 
biomarkers  is  weak  and  suggests  that  no  single  biomarker  can 
fully  explain  the  behavior  of  another.11  W  hen  adjusting  for  the 
effect  of  time,  the  vessel  diameter  was  found  to  be  associated 
with  the  redox  ratio  of  the  adjacent  tissue,  which  indicates 
that  vessel  diameter  plays  a  key  role  in  tumor  metabolic  ac¬ 
tivity  and  cycling  hypoxia  in  tumors.43  Future  studies  will 
exploit  these  biomarkers  to  investigate  functional  changes  in 


Fig.  5  Normalized  redox  ratio  as  a  function  of  distance  from  the  near¬ 
est  vessel  for  tumors  (compiled  from  n  =  94  line  profiles  corresponding 
to  the  94  regions  of  interest  from  the  three  redox  animals  across  the 
entire  tumor  time  course),  and  for  normal  tissues  (compiled  from  n 
=  36  line  profiles  from  six  regions  of  interest  from  each  of  six  normal 
animals).  All  profiles  were  normalized  to  the  point  closest  to  the  ves¬ 
sel  before  averaging  across  all  line  profiles.  For  the  tumors,  there  is  a 
significant  (p<0.05)  decrease  in  the  redox  ratio  between  the  point 
closest  to  the  vessel  (point  zero)  and  all  points  greater  than  30  pm 
from  point  zero.  Normal  tissues  did  not  show  a  significant  increase  or 
decrease  in  the  redox  ratio  as  a  function  of  the  distance  from  the 
nearest  vessel.  Error  bars  are  standard  error. 


oxygen  supply  and  demand  in  response  to  therapeutic  strate¬ 
gies.  For  example,  it  is  possible  that  antiangiogenic  therapies 
may  stabilize  cycling  hypoxia  by  reducing  the  magnitude  of 
changes  in  vascular  diameter  and  blood  flow  rate  over  time.3 
It  is  anticipated,  however,  that  antiangiogenic  therapies  would 
have  less  influence  on  redox  ratio.  The  combination  of  tech¬ 
nologies  presented  here  would  permit  a  direct  test  of  this  hy¬ 
pothesis. 

Extracted  biomarkers  including  the  hemoglobin 
saturation,44,45  maximum  velocity  and  vessel  diameter,14  and 
shear  rate,46  are  all  in  qualitative  agreement  with  previous 
studies  that  employed  different  techniques  in  window  cham¬ 
ber  tumors,  as  well  as  a  previous  study  that  employed  the 
same  techniques  in  a  different  window  chamber  tumor 
model.28  Variations  in  red  blood  cell  flux  and  p02,15  hemo¬ 
globin  oxygen  saturation,23  and  diameter44  have  previously 


Table  1  Pearson  correlation  coefficients  for  hemoglobin  saturation 
(FHbSat)  versus  maximum  velocity  (V^),  flow,  diameter,  and  shear 
rate  (top  line)  and  for  the  normalized  redox  ratio,  normalized  hemo¬ 
globin  saturation,  normalized  maximum  velocity,  normalized  flow, 
normalized  diameter  and  normalized  shear  rate  (bottom  line).  The 
data  are  compiled  from  all  regions  of  interest  and  time  points  grouped 
together,  with  n=143  observations  for  the  top  line  and  n  =  94  obser¬ 
vations  for  the  bottom  line.  The  correlation  coefficient  is  listed  for 
comparisons  with  p<0.05. 


v 

r  max 

(mm/s) 

Flow 
(mm3/ s) 

Diameter 

(pm) 

Shear  rate 

(1/s) 

HbSat  (%) 

0.49 

0.37 

p>0.05 

0.49 

HbSat 

y  max 

Flow 

Diameter 

Shear  rate 

Redox 

0.50 

0.28 

0.32 

0.43 

p>0.05 
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Fig.  6  Longitudinal  variations  in  all  biomarkers  as  a  function  of  time. 
Values  at  each  time  point  were  averaged  across  all  regions  of  interest 
and  animals  for  the  hemoglobin  saturation  (HbSat),  flow,  diameter, 
maximum  velocity  (Vmax),  and  shear  rate  (n= 21),  as  well  as  for  the 
redox  ratio  (n=  14).  Note  that  the  36-h  time  point  has  four  fewer 
samples  (see  Sec.  2).  The  linear  mixed  model  determined  that  all 
biomarkers  significantly  change  (p<  0.05)  with  time,  and  that  the  ves¬ 
sel  diameter  is  associated  with  the  redox  ratio  of  the  adjacent  tissue 
(p<0.05)  after  adjusting  for  the  effect  of  time.  Error  bars  are  standard 
error. 


been  observed  in  window  chamber  tumors  over  time  scales 
ranging  from  one  hour  to  several  days. 

Previous  studies  have  employed  different  methods  to  in¬ 
vestigate  the  relationship  between  some  of  the  biomarkers 
presented  in  this  work.  Intravascular p02  and  blood  flow  rate, 
as  well  as  interstitial  p02  and  red  blood  cell  velocity  or  vessel 
diameter,  were  found  to  have  no  correlation  in  a  previous 
study.13  Helmlinger  et  al.  used  phosphorescence  lifetime  im¬ 
aging  to  quantify  p02  and  video  microscopy  to  quantify  red 
blood  cell  velocity  and  vessel  diameter.  Although  hemoglobin 
oxygen  saturation  and  p02  are  related,  they  are  not  the  same 
endpoint,  and  video  microscopy  provides  a  2-D  projection 
rather  than  an  absolute  measure  of  blood  velocity.  The  study 
design  and  tumor  models  employed  in  the  previous  and  cur¬ 
rent  studies  also  differ.  Previous  studies  employing  the  same 
techniques  as  the  current  study  have  found  a  significant  cor¬ 
relation  between  the  hemoglobin  saturation  and  blood  veloc¬ 
ity,  blood  flow,  and  shear  rate  in  a  different  tumor  model.28 
Additional  studies15  indicate  that  vessel  remodeling  with  tu¬ 
mor  growth  alters  flow  distributions,  and  the  relationship  be¬ 
tween  blood  flow  and  tissue  oxygenation  in  tumors  is  hetero¬ 
geneous. 

Demand  for  oxygen  in  tumor  cells  is  created  by  the  normal 
respirative  process  of  the  cells.32  The  redox  ratio  is  a  measure 
of  cell  respiration  (oxygen  demand),30  while  the  hemoglobin 
oxygen  saturation  is  one  component  of  oxygen  supply.  There¬ 
fore,  it  is  intuitive  that  the  hemoglobin  saturation  of  a  vessel  is 
somewhat  related  to  the  redox  ratio  of  the  surrounding  tissue. 
However,  we  do  not  expect  a  one-to-one  correlation  between 
the  redox  ratio  and  hemoglobin  saturation.  Oxygen  delivery  is 
affected  by  the  hemoglobin  saturation  as  well  as  the  vascular 
architecture,  hematocrit,  and  blood  flow.32  The  redox  ratio  is 
indicative  of  changes  in  oxygen  supply  as  well  as  metabolite 
availability  and  mitochondrial  membrane  potential.47  Thus, 
the  redox  ratio,  blood  flow  and  hemoglobin  saturation  provide 
complementary  information. 

Figure  5  shows  a  decrease  in  the  redox  ratio  with  increas¬ 
ing  distance  from  the  nearest  vessel  in  tumors,  and  no  such 
decrease  in  the  redox  ratio  with  distance  from  the  nearest 
vessel  in  normal  tissues.  We  hypothesize  that  this  tumor  gra¬ 
dient  is  due  to  increased  metabolic  activity  coupled  with  com¬ 


promised  microvascular  oxygen  transport  in  tumors  compared 
to  normal  tissues.  Future  studies  will  investigate  this  hypoth¬ 
esis. 

In  this  study,  absorption  and  fluorescence  microscopy  were 
combined  with  OCT  to  unravel  the  contributions  of  blood 
flow,  blood  oxygenation,  and  metabolic  demand  on  cycling 
hypoxia.  These  findings  could  guide  the  design  and  develop¬ 
ment  of  therapeutic  strategies  to  mitigate  cycling  hypoxia,  and 
thus  reduce  tumor  resistance  to  radiation  and  chemotherapies. 
Second,  this  study  experimentally  demonstrates  the  attributes 
of  our  multifunctional  optical  imaging  approach  for  mapping 
oxygen  supply  and  demand  in  vivo.  This  approach  has  broad 
applicability  in  the  study  of  other  diseases  that  are  affected  by 
oxygen  supply  and  demand,  such  as  cerebral  hypoxia  (stroke), 
cardiac  ischemia,  and  Alzheimer's  disease.  Moreover,  these 
processes  can  be  studied  with  multifunctional  optical  imaging 
in  vivo  and  do  not  require  contrast-enhancing  dyes. 
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